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Abstract: We show that Af = 1 gauge theories with an adjoint chiral multiplet 
admit a wide class of large- A^ double-scaling limits where A^ is taken to infinity in 
a way coordinated with a tuning of the bare superpotential. The tuning is such 
' that the theory is near an Argyres-Douglas-type singularity where a set of non-local 

dibaryons becomes massless in conjunction with a set of confining strings becoming 
^ ■ tensionless. The doubly-scaled theory consists of two decoupled sectors, one whose 

' spectrum and interactions follow the usual large- A^ scaling whilst the other has light 

states of fixed mass in the large- A^ limit which subvert the usual large- A^ scaling and 
lead to an interacting theory in the limit. F-term properties of this interacting sector 
can be calculated using a Dijkgraaf-Vafa matrix model and in this context the double- 
scaling limit is precisely the kind investigated in the "old matrix model" to describe 
two-dimensional gravity coupled to c < 1 conformal field theories. In particular, the 
old matrix model double-scaling limit describes a sector of a gauge theory with a 
mass gap and light meson-like composite states, the approximate Goldstone boson 
of super conformal invariance, with a mass which is fixed in the double-scaling limit. 
Consequently, the gravitational F-terms in these cases satisfy the string equation of 
the KdV hierarchy. 
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1. Introduction 



't Hooft argued many years ago [1] that the perturbative expansion of theories in- 
volving matrices has an interesting interpretation when N, the size of the matrices, 
becomes large. The idea is to first re-write the perturbative expansion in Qg, the 
coupling, and N, in terms of 5* = QsN, the '"t Hooft Coupling", and Qg- The Feyn- 
man graphs of the theory can be sorted on the basis of the topology of an imaginary 
two-dimensional surface on which the diagrams can be drawn. In terms of {S,gs), 
the explicit Qs dependence of a graph is determined solely by the topology of the 
imaginary surface: if g is the genus then the dependence is gl^~^. At a given genus, 
there is then a perturbative series in the 't Hooft couphng S. 

An example of such matrix theories are Yang-Mills theories. These ideas led 
't Hooft to suggest that Yang-Mills theories may have an alternative description 
in which the imaginary two-dimensional surfaces become physically realized as the 
world-sheets of a string theory. In this picture, Qs = Qym becomes the string coupling 
weighting each handle in string world-sheet perturbation theory, so that the string 
theory is weakly coupled when N is large. Quite how gauge theories are equivalent 
to a theory of strings has taken many years to understand in a concrete way. 

The topological expansion of matrix field theories is quite general and applies, 
for example, to simple matrix theories (zero-dimensional field theories), where there 
is no spacetime. The fact that such simple theories can give rise to two-dimensional 
surfaces has been an idea of great import. In the "old matrix model" epoch the idea 
was to use matrix models to describe models of two-dimensional gravity coupled 
to various matter systems. Here, once again the surfaces arise in an auxiliary way 
from the expansion of the matrix model in the string coupling. In these theories the 
way that the Feynman graphs give rise to a two-dimensional surface is completely 
transparent. The 't Hooft expansion of, say, the free energy has the form 

oo 

F^J2^^(09'r' . (1.1) 

9=0 

Here, ^ label a general set of parameters in the theory (suitably scaled in N), in- 
cluding the 't Hooft couphng. The surfaces arise when one takes the large- A?" limit 
with fixed 't Hooft coupling, as before, but, in addition, one takes a parameter ^ to 
a critical value. At the critical point ^ = ^o, the contributions from a given genus 
Fg{^) diverge in a particular way: 

- \C - Col"^'-''^ , (1.2) 

where 77 is a critical exponent which characterizes the universality of the critical point. 
Intuitively, near the critical point the Feynman graphs of a given genus become very 
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dense and describe a surface. The coordinated, or double-scaling limit, involves 
taking [2] 

N^oo, C^Co, A = 7V|e-eor = fixed. (1.3) 
In this limit there is an effective string coupling 

{g,U ~ (1.4) 
which weights the contributions from different genera. 

What makes the old matrix model story so compelling is the ease at which the 
two-dimensional surfaces arise. This is not so clear in Yang-Mills theories where the 
surfaces are thought to be associated to a dual string description. Realizing this 
picture explicitly was brilliantly achieved by Maldacena in the case of maximally 
supersymmetric J\f — 4 gauge theories [3]. In this case the dual string theory is 
a ten-dimensional Type IIB superstring on the space AdS^ x with background 
Ramond-Ramond (RR) flux. Unfortunately such string theories have proved very 
difficult to solve because of the presence of the RR flux. It would be very useful to 
have a situation where the string theory appeared in a very simple way as in the 
old matrix model. What we have in mind is a realization of a double-scaling limit 
in a Yang-Mills theory where is taken to infinity as certain couplings in a bare 
potential are taken to critical values in such a way that the string dual naturally 
arises. Double-scahng hmits of this kind have been recently studied in [4] for a class 
of jV = 1 gauge theories in a partially confining phase. There it was proposed that 
they admit a dual description in terms of a non-critical string theory of the type first 
introduced in [5]. The worldsheet theory corresponding to this string background 
is exactly solvable and there are no RR fiuxcs. The gauge theories involve a gauge 
group G = U {N) which is confined down to an abelian group. The theories have an 
M = 1 vector multiplet and an adjoint- valued chiral multiplet $ (precisely the fields 
of an A/" = 2 theory).^ There is an explicit bare superpotential 

t+i 

H/bare(*) = Tr W {^) , W {x) = nY^x' . (1.5) 

The prefactor ensures conventional large- scaling. Classically, each of the eigenval- 
ues of the scalar component of $ can be one of the i critical points of W{x), 

£ 

W'{x) = Ne l[{x -ai) , e = ge+i . (1.6) 

i=l 

^The /3-deformation of AA = 4 Super Yang-Mills also has a partially confining phase that has 
been studied in [6-8] . A double-scaling limit and the relative non-critical string dual have also been 
proposed in [4]. 
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We define the multiplicity of eigenvalues at the critical point a, as Ni > 0. At least 
in an appropriate limit, this leads to a classical Higgs effect U{N) Y[l=iUi^i)^ 
where Yli=i ~ ^- Each U{Ni) will then subsequently confine to leave a U{1) 
factor in the JR. We will denote by s the number of A^j > 0, so that the IR gauge 
group is Uiiy . Of course, such a picture of a classical Higgs effect plus confinement 
is not expected to be valid for all the values of the parameters but is a useful limit to 
have in mind. In particular, vacua with the same rank of IR group, but with different 
"filling fractions" Ni/N, are known, under certain circumstances, to be continuously 
connected in the space of superpotentials [9] [10]. 

The spectrum of these partially confining theories is expected to be rather compli- 
cated and since there is only Af = 1 supersymmetry and no notion of a BPS state, we 
have little analytic control over particle masses. However, we expect that for a generic 
vacuum the conventional large- picture applies (see for example [11]): as well as 
the massless abelian sector, there will be a tower of glueballs and " dimesons" , formed 
by one "diquark" Q^s transforming in the (A,., A^,) representation of U{Nr) x U{Ns) 
and the corresponding anti-diquark Qrs, whose masses are A^ independent. There 
will also be dibaryons, whose masses scale like A^. The interactions of the glueballs 
and mesons are suppressed by l/N,"^ whereas the interactions between the baryons 
and mesons are not suppressed. However, if there is a point i^o in the space of pa- 
rameters {gi} where some state becomes massless then the usual large- A^ picture can 
break down. These states are baryonic (or composites of these) with a mass that 
generically goes like A^, but near the critical point one finds that M ~ A^|^ — 
for some exponent t]. In the double-scaling limit in which N ^ oo and ^ — > 
such that M is fixed, the theory decomposes into two sectors H 7Y_ x T-C+. For 
the sector 7i„, the usual rules of the large- A^ expansion are modified by diagrams 
involving the state of mass M running around in loops. As described in [4] there 



is now an effective string coupling as in (|1.4| ). The sector on the other hand 



decouples and is bound by conventional large- A^ reasoning (and hence is not inter- 
esting). What is remarkable is that there is a natural candidate for the string dual 
to the sector 7i_, namely a non-critical string background of the type introduced 
in [12]. These non-critical backgrounds were originally studied as holographic du- 
als to certain Little String Theories (LST) and double-scaled LST [13-16]. In the 
low-energy limit these LST reduce to 4d A/" = 2 theories in the proximity of Argyres- 
Douglas singularities which are non-trivial superconformal field theories [17-19]. The 
fact that these backgrounds preserve M = 2 supersymmetry seems to contradict the 



^The interactions between mesons are suppressed by integer powers of 1/iV instead of the usual 
because these states are formed by one "diquark" Q^s transforming in the {Nr,Ns) repre- 
sentation of U{Nr) X U{Ns) and the corresponding anti-diquark Qrs- Therefore the three-point 
couphng of these "dimesons" gains a suppression factor ~ 1/VN from each of the two groups, 
Nr^N.^N [4]. 
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proposed duality. Remarkably, the F-terms of the theory are consistent with this 
supersymmetry enhancement in the 7i_ sector [4]. It is interesting to try to extend 
this kind of duality and search for more general classes of double-scaling limits in 
supersymmetric gauge theories. In this paper, we will achieve this aim and we will 
show in many cases that these double- scaling limits exist and that the Hilbert space 
of the theory still exhibits the above splitting into two decoupled sectors. As be- 
fore, the sector has non-trivial dynamics in the double- scaling limit weighted by 
the effective string coupling (|1.4| ). However, in the generic situation we consider, 
this sector only has M = 1 supersymmetry since the effects of the superpotential 
do not vanish in the limit. This is different from the case studied in [4], where the 
superpotential vanished in the Ti- sector, which is consistent with the enhancement 
to A/" = 2 supersymmetry. Consequently, what is missing from our analysis is the 
precise identification of the dual string theory itself ^. We will also find examples 
where the sector 7Y_ does not contain any U{1) group but has a mass gap and light 
meson-like composite states. 



So our focus is on large- iV double-scaling limits in supersymmetric gauge theories. 
As we have argued, in order to have such a limit we need a scenario in which there 
are states which we can tune to have a finite mass as N ^ oo. The question is 
how can we verify that such double-scaling limits actually exist because we do not 
generally have analytic control over the masses of particles in an A/" = 1 theory? 
The answer is that, under certain circumstances, we can see the effects of a nearly 
massless particle in various F-terms and these can often be calculated exactly. For 
instance, if the particle is charged under a f/(l) gauge symmetry that survives in the 
IR, then we will see the effects of the light particle in the associated renormalization 
of the U{1) coupling. This was one of the strategies employed in [4] to explore these 
double-scaling limits. In the present work, we will find situations where the light 
states are not charged under any U{1) in the IR. So in order to "see" these kinds 
of states we will have to investigate other F-term couplings of the gauge theory. 
The idea is to consider such couplings that arise when the gauge theory is coupled 
to supergravity. In that case there is a whole set of gravitational F-terms [21,22]. 
There are two series of terms that interest us 



and 




construction of four-dimcnsional A/" = 1 non-critical strings dual to A/" = 1 quiver gauge 
theories has recently appeared in [20] 
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Here, Wais-y and Fap are the Af = 2 gravitino and gravi-photon superfields of the 
gravity sector, respectively. The Fg{S) are the gauge theory objects: functions of 
the s glueball superfields {Si}. Note that the g = term of Fi is the ordinary 
glueball superpotential: 

i=l * 

which determines the vacua of the theory.^ 

In the usual large- limit, the contribution from a genus g graph is suppressed by 
j\j2~2g ^ However, if there is a scenario as described above where a state is becoming 
massless then this suppression can be subverted. The effect of such a state on Fg 
can be crudely estimated by considering the light state propagating around a 1-loop 
graph with g insertions [4,24]. This suggests Fg ~ M^"^^, which together with the 
identification M ~ A means that the gravitational F-term Fg behaves just like the 
genus g free energy of a matrix model as in (|1.2|) . The conclusion is that the quantities 
Fg appearing in the gravitational F-terms provide another signal of a double-scaling 
limit. Of course it is no accident that these couplings behave like the free energy of a 
matrix model since they can be calculated by just such a matrix model as described 
by Dijkgraaf and Vafa [25]. This matrix model is defined in terms of an iV x iV 
matrix $ (A^ here is not to be confused with of the field theory) with the bare 
superpotential in ( |L5D : 

oo „ 

exp ^ / exp (- (7;iTrH^(<l>)) . (1.10) 

The integral is to be understood as a saddle-point expansion around a critical point 
where iVj of the eigenvalues sit in the critical point a, (but only if the associated 
Ni > 0) and by definition the glueball superfields are identified with the quantities 

s 

S, = gA, S = J2S^ = 9sN (1.11) 

i=l 

in the matrix model. It is worth remarking that in the matrix model we work at 
large and that this procedure would yield exact results for any value of A^. But, 
of course, we also want to take A^ to be large in the double-scaling limit. 

It should now be apparent that the double-scaling limit of the gravitational F- 
terms ( |1.7| ) ( |1.8| ) in the gauge theory corresponds precisely to the notion of a double- 
scaling limit in the "old" matrix model and in certain cases we can use standard 
matrix model technology to investigate the phenomenon. 

^We are assuming that the VEVs of the supergravity fields vanish otherwise we have to keep all 
the terms in Fi in the superpotential as in [23]. 
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The results presented in this paper have some overlap with earlier work. The 
idea that the standard large- expansion breaks down near the critical points in the 
moduli/parameter space of four-dimensional gauge theories has been emphasized in 
a series of papers by Ferrari [26-28] , which also contain proposals for double-scaling 
limits. The four-dimensional interpretation of double-scaling in the Dijkgraaf-Vafa 
matrix model has also been discussed before in [25,26]. 

The plan of the paper is as follows. In section ^, we will give a general description 
of double-scaling limits and explain how they can be realized on-shell. In section |], 
we will study the exact F-term couplings of the various models and deduce that in 
the double-scaling limit the Hilbert space of theory decomposes into two decoupled 
sectors 7i+ and where Ti- keeps finite interactions in the limit. In section ^, we 
will analyze the behaviour of the higher genus terms of the matrix model free energy 
Fg at the critical point and show that Fg ~ A^~^^. In section ^, we will discuss in 
detail the case where the near-critical curve has double points. We will argue that in 
the critical limit where a double point collides with a branch point a meson-like state, 
which can be thought of as the bound state of two dibaryons, becomes massless. 



2. Engineering Double-Scaling Limits in the Matrix Model 



The central object in matrix model theory is the resolvent 

^(x) = iTr^-^. (2.1) 

At leading order in the 1/A^ expansion, lj{x) is valued on the spectral curve S, a 
hyper-elliptic Riemann surface 

The numerical pref actor is chosen for convenience. In terms of this curve 

ujix) = ^{W'ix)-Neyix)) . (2.3) 

In ( p.2| ), fe-i{x) is a polynomial of order i—1 whose i coefficients are moduli that are 
determined by the Si. In general, the spectral curve can be viewed as a double-cover 
of the complex plane connected by i cuts. For the saddle-point of interest only s of 
the cuts are opened and so only s of the moduli /^_i(x) can vary. Consequently y{x) 
has 2s branch points and i — s zeros: ^ 

S : y^ = Z„(x)V2,(x) (2.4) 



^Occasionally, for clarity, we indicate the order of a polynomial by a subscript. 
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where £ = m + s and 

m 2s 

Zm{x) = Y\_ix - Zj) , a2s{x) = Y\_ix - aj) . (2.5) 

The remaining moduh are related to the s parameters {Si} by ( |1 . 1 ![ ) 

S, = gsNi = Ne (f y dx , (2.6) 

J A, 

where the cycle Ai encircles the cut which opens out around the critical point of 
W{x). 

Experience with the old matrix model teaches us that double-scaling limits can 
exist when the parameters in the potential are varied in such a way that combinations 
of branch and double points come together. In the neighbourhood of such a critical 
point, ^ 

n 

CZ^{xfBn{x) , 5„(x) = l[{x - h) , (2.7) 

i=l 

where Zj^hi —>■ xq, which we can take, without loss of generality, to be Xq = 0. The 
double-scaling limit involves first taking a ^ 

X = ax , Zi = azi , bj = abj (2.8) 

while keeping tilded quantities fixed. In the limit, we can define the near-critical 
curve T,_-J 

S_ : yl = Z^xYBr^ix) . (2.9) 

We will argue, generalizing [4], that in the limit a — > 0, in its sense as a complex 
manifold, the curve S factorizes as S_ U S+. The complement to the near-critical 
curve is of the form 

S+ : yl = x'"'+^C2s-n{x) . (2.10) 
where C2s-n{x) is regular at a = 0. 

In the a — > limit, we will show in Section ^ that the genus g free energy gets a 
dominant contribution from S_ of the form 

Fg ~ (iVa(™+"/2+i))'-'^ . (2.11) 

Note that in this equation N is the one from the field theory and not the matrix 
model N. This motivates us to define the double-scaling limit 

a^O, iV^oo, A = Ara"^+"/2+i = const . (2.12) 

^We have chosen for convenience to take all the double zeros {zj} into the critical region. 
''For polynomials, we use the notation f{x) = Yiii^ ~ fi)^ where f{x) = Yiii^ ~ fi)i x = ax and 
fi = afi. 
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Moreover, the most singular terms in a in ( |2.11 ) depend only on the near-critical 
curve (|2.9| ) in a universal way. 



2.1 Engineering the double-scaling limit on-shell 

However, there is still an important issue to address. In the context of supersym- 
metric gauge theories, the moduli {Si} are fixed by extremizing the glueball su- 
perpotential ( |1.9|) . It is not, a priori, clear whether a double-scaling limit can be 
reached whilst simultaneously being on-shell with respect to the glueball superpo- 
tential. We now address this issue and show that suitable choices of the coupling 
constants {gi\ do indeed allow for a double-scaling limit on-shell with respect to the 
glueball superpotential. In general, the potentials required are non-minimal. How- 
ever, this is irrelevant for extracting the universal behaviour that only depends on 
the near-critical curve (|2.71 ). 

So the problem before us is to show that the critical point can be reached simul- 
taneously with being at a critical point of the glueball superpotential. It is rather 
difficult to find the critical points of the latter directly. Fortunately another more 
tractable method consists of comparing the matrix model spectral curve (|2.2|) , the 
"TV = 1 curve" , with the Seiberg- Witten curve of the underlying M = 2 theory that 
results when the potential vanishes. The latter has the form 

yl^ = PM{xf - AK''' , (2.13) 

where Pn{x) = YliLii^ ~ 4>i)- Here, are a set of coordinates on the Coulomb 
branch of the M = 2 theory and A is the usual scale of strong-coupling effects in the 
M = 2 theory. 

When the Af = 2 theory is deformed by addition of the superpotential (|1.5|), it 
can be shown that a vacuum exists when the Seiberg- Witten curve and the A/" = 1 
curve represent the same underlying Riemann surface [10,29,30]. In concrete terms 
this means that, on-shell, 

y' = {Wixf + fs^m^iix)) = ZUxfaUx) , ^^-^^^ 

In these equations, H]\f_s{x), a2s{x), Zm{x) are polynomials of the indicated order, 
and we choose (in order to remove some redundancies) 

Hn.s{x) = x^-^ + ■ ■ ■ , a^six) = + ■ ■ ■ , Zrn{x) = + ■ ■ ■ . (2.15) 
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Both curves describe the same underlying Riemann surface, namely the reduced 
curve of genus s — 1 which is a hyper-elliptic double-cover of the complex plane with 
s cuts. All-in-all there are 2(A^ + /) equations for the same number of unknowns 
in {P, H,a, Z, f}. There are many solutions to these equations and we can make 
contact with the description of the vacua in Section 1 by taking the classical limit 
A 0; whence 

e e 
Pn{x) ^ - ^*)'^' ' 5Z AT, = AT , (2.16) 

i=l i=l 

SO Aj of the eigenvalues of the Higgs field classically lie at the critical point of 
W{x). Quantum effects then have the effect of opening the points into cuts (if 
A'i > 0). The number of Ni > 0, i.e the number of cuts, is equal to s = i — m. 

We now to turn to explicit solutions of ( |2.14|) . The method we shall adopt is 
to first find solutions for a U (p) gauge theory and then apply the "multiplication by 
N/p map" [30], with N/p integer. This will yield a solution for a U{N) gauge group 
and will allow to take a large A^ limit with p fixed. 

2.2 No double points 

We now describe how to engineer the case where the near critical curve ( p.7|) has no 
double points, so m = 0. This is the situation considered in [4,31,32]. In this case, 
we first consider the consistency conditions (|2.14]) for a U{p = n) gauge theory with 
W{x) of order £ = n + 1. In this case, ( p.l4|) are trivially satisfied with 

W'{x) = NePnix) , /„_i(x) = -ANh^A^^ . (2.17) 

Notice that with our minimal choice of potential, the on-shell curve actually implies 
that S = since the coefficient of x""^ in fn~i{x) vanishes and so the resolvent falls 
faster than 1/x at infinity. This, of course, is pathological from the point-of-view of 
the old matrix model and may be remedied by using a non-minimal potential with 
extra branch points or double points outside the critical region. However, in the 
holomorphic context in which we are working, having S" = is perfectly acceptable 
and we stick with it. The on-shell curve consists of an n-cut hyperelliptic curve and 
one can verify, by taking the classical limit, that = 1, i = 1, . . . ,n. The double- 
scaling limit involves a situation where n branch points, one from each of the cuts, 
come together. This can be arranged by having 

n 

W'{x) = Ne{Bn{x) + 2A") , Bn{x) = \{{x - bj) (2.18) 
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and then taking the hmit (2^). In this case, the near-critical curve S_ (2^) is of 
the form 

yl = Bn{S:) . (2.19) 

The important point is that we can tune to the critical region whilst keeping the 
theory on-shell with respect to the glueball superpotential by simply changing the 
parameters which appear in the potential. 

Now that we have found a suitable vacuum of a U{n) theory, we now lift this to 
a U (N) theory with the multiplication by N/n map [30]. Under this map, the A/" = 1 
curve remains intact, including the potential W{x) whilst the Seiberg-Witten curve 
of the U{N) theory is 

ylw = PN{xf - 4A^^ = A^(^-")W^_^ (^)'(^"^^^' ~ ^^'"^ ' ^^-^^^ 
where Un_ {x) is a Chebishev polynomial of the second kind. The vacuum of the 

n 

U (N) theory has Ni = N/ n, i = 1, . . . ,n. 

Notice that in the near critical region the Seiberg-Witten curve is identical to 
E_, up to a rescaling: 

ylw^(—y^'''-"Bn{x). (2.21) 



v n 



This is simply a reflection of the observation of [4] that the decoupled sector has 
enhanced M = 2 supersymmetry. Moreover, if C is a cycle which is vanishing as 
a — > then the integral of the Seiberg-Witten differential around C, which gives the 
mass of a BPS state carrying electric and magnetic charges in the theory, becomes 

r xP^ix) dx _^ _^_„/2^^„/2+i ly_di. (2.22) 
Jc Vsw Jc 

Notice that in the double- scaling limit ( p.l2|) (with m = 0) the mass of the state 
is fixed. This state is a dibaryon that carries electric and magnetic charges of the 
IR gauge group. In the double-scaling limit, therefore, a set of mutually non-local 
dibaryons become very light. ^ In fact, the Seiberg-Witten curve at the critical point, 
a = 0, has the form 



?/sw = 4(-) A^^--x^ (2.23) 
which describes a Z„ or Argyres-Douglas singularity [17-19]. 



^For n — 2 there is only a single light dibaryon. 
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2.3 With double points 



For the case with double points, we cannot simply take two of the branch points 



{bj} in (|2.18| ) above to be the same. If we simply did that then the zero of the 



Seiberg-Witten curve, by which we mean a zero of the polynomial H^-s in ( |2.14| ) 



would also be a zero of the J\f = 1 curve as well. By the analysis of [29], this would 
imply that the condensate of the associated massless dibaryon would be vanishing. 
On the contrary, if the zero of the Seiberg-Witten curve were not a zero of the A/" = 1 
curve, the putative massless dibaryon would be condensed and the dual U{1) would 
be confined. We need to arrange the situation so that any zero of the Seiberg-Witten 
curve is not simultaneously a zero of the A/" = 1 curve. 

A suitable Af = 1 curve which reduces to (|2.7|) in the near-critical region is 

= ZUxfB^ix) {B^{x)H,{xf + AK^'^^) . (2.24) 

In this case, we have i = m + n + r, s = n + r and 

W'{x) = NeZ^{x)B^{x)Hr{x) , /,_i(x) = AN^e^A^'+^Z^xfE^ix) . (2.25) 

Notice that in order that fe-i{x) has order less than i we require r > m. The curve 
( p.24| ) is actually on-shell with respect to the Seiberg-Witten curve of a f/ (2r + n) 
theory with 

P2r+n{x) = Hr{xfBn{x) + 2A2'-+" . (2.26) 

In the classical limit, we have two eigenvalues at each of the zeros of Hj.{x) and one 
in each of the zeros of Bn{x). Once again we can employ the multiplication map 
( |2.20|) (with n replaced by 2r + n) to find the vacuum of the U{N) theory we are 
after. 

Notice that the double points of the Seiberg-Witten curve {hi} are not generally 
zeros of the curve ( ^.24[ ), which means that the associated dyons are condensed. The 



near-critical curve S_ in this case is 

= Zm{ifBn{i) , (2.27) 
while in the near-critical region the Seiberg-Witten curve becomes 

^ i7,,(0)2i?„(x) . (2.28) 

where we assumed that the zeros of Hr{x) lie outside the critical region. In this case, 
the integral of the Seiberg-Witten differential around a vanishing cycle diverges in 
the double-scaling limit: 

"'^Ni'')^'' ^ Nan/2+1 ^ ^^-m ^ ^ _ ^2.29) 

c ysw 
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So in contrast to the case with no double points, the dibaryon states are very heavy. 
In addition, the dyon condensate associated to the zero hi of Hr{x) is given by an 
exact formula [29] 

{niirhi) = Ney{hi) ^ N ^ oo , (2.30) 

where we have assumed that hi stays fixed as a ^ 0. So in the double- scaling limit 
the value of the condensate and hence the confinement scale in the dual U{1), or 
string tension, occurs at a very high mass scale. 

We now have a puzzle. How can there be an interesting double-scaling limit 
in the gauge theory in this case if there are no light dibaryon as in the previous 
example? The answer is that there are other light mesonic states in the theory with 
a mass ~ A that we will identify in Section ^. 

Notice that contrary to our choice above, if we scale /ij ^ as a — then the 
tensions of the confining strings vanish and the theory is at an A/" = 1 superconformal 
fixed point in the infra-red corresponding to one of the A/" = 1 Argyres-Douglas-type 
singularities described in [31]. As the double points of the Seiberg-Witten curve hi 
move away from the associated dyons condense and the superconformal invariance 
is spontaneously broken. The resulting Goldstone modes will play an important role 
in the ensuing story and resolve the puzzle alluded to above. 



3. The Double- Scaling Limit of F- Terms 



Before we consider the a — limit of the free energy, it is useful to the consider 
this limit for other F-terms in the low-energy effective action. The effective action is 
written in terms of chiral superfields Si and Wai which are defined as gauge-invariant 
single-trace operators [33] 

Wnl = (t dx 

2m Air 

It will also be convenient to define component fields for each of these superfields, 

Si = si + e^x? + ■ ■ ■ , w^i = X^i + Opf^i + ■■■ . (3.2) 

The component fields, si and /; are bosonic single trace operators whilst xi ^"^^ 
are fermionic single trace operators. In the large- A^ limit, these operators should 
create bosonic and fermionic colour- singlet single particle states respectively. It is 



X - $ 



X - $ 
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instructive to consider the interaction vertices for these fields contained in the F-term 
effective action whose general form is given by [25, 34, 35] 



£^ = Im / £d (lygb + W^eff ) 



(3.3) 



where 



k,l 



(3.4) 



Expanding (p.3|) in components on-shell, we find terms like 



(3.5) 



where 








(3.6) 



ds,,ds,,...ds. 



for L = 2, 3, 4. In the large- limit, V^^'^ scales like iV^ ^ . We will also consider the 
2-point vertex coming from the glueball superpotential 



The matrix Hi- therefore effectively determines the masses of the chiral multiplets 
Si. Note that, in the large- hmit, H^'^^ scales like iV°. 

We begin by considering the couplings V^J''^ of the low-energy f/(l)* gauge group. 
Each of the ?7(l)'s is associated to one of the glueball fields Si, or equivalently the 
set of 1-cycles {Ai\ on S. If we ignore the f/(l) associated to the overall 't Hooft 
coupling S, or the cycle A^o = Tl,l=i "which can be pulled off to infinity, the 
couplings of the remaining ones are simply the elements of the period matrix of S. 
In order to take the a — >■ limit, it is useful to choose a new basis of 1-cycles {Aj, 
i = 1, . . . , s — 1, which is specifically adapted to the factorization S S_ U S+. 
The subset of cycles with i = 1, . . . ,[n/2] vanish at the critical point while the cycles 
i = [r;,/2] + 1, . . . , s — 1 are the remaining cycles which have zero intersection with all 
the vanishing cycles. 

If we define the periods on S 




(3.7) 



where 




dS.dSj ■ 



(3.8) 




(3.9) 
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then the period matrix, in this basis, is simply 

U = N~^M. (3.10) 

In the Appendix, we calculate the a — » limit of these matrices. The results are 
summarized in ( |A.4| ) and ( |A.7|) . Using these results, we have 




Nz'_M^^l+AfM^^l\ 



(3.11) 



Let us look more closely at the structure of each block in the above matrix. First of 
all, by (Q 

(^--)., ~ if > (^-)., - /f ^(^0 , (3.12) 

which implies 

{N^^):^ (M__)^.^ = fip--\k)f^f{bi) = . (3.13) 

Furthermore, since Nz}_ vanishes in the limit a — > 0, we find that 

Ar-lMi°j + ATmI"! = NZ'_ M^'l - iVrl Ari°j(Ar|°j)-'M|°j ^ . (3.14) 

Therefore, the period matrix has the following block-diagonal form in the double- 
scaling limit 

The upper block 11^ is actually the period matrix of the near-critical spectral curve 
S_ (|2^ ) since the cycles {Ai,Bi}, for i < [n/2] form a standard homology basis for 
S„. Similarly, the lower block 11^ is the period matrix of S+. So in the limit a —>■ 
the curve S factorizes as S_ U S+. The fact that the period matrix factorizes is 
evidence of the more stringent claim that the whole theory consists of two decoupled 
sectors and 7i+ in the double-scaling limit. Note that although we did not 
consider it, the U{1) associated to S only couples to the 7i+ sector. 

We can extend this discussion to include other F-terms that are derived from 
the glueball superpotential. For example, consider the 3-point vertex 

VS! = (3.16) 



dS^dSjdS, 



k 



Here, the Si as defined as in ( |1.11D but with respect to the cycles Ai. They are 



related to the Si by an electro-magnetic duality transformation. There is a closed 
expression for these couplings of the form [36-38] 

^.S'-^Ef^e.,™. (3.17) 
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where {ooj} are the holomorphic 1-forms normahzed with respect to the basis {Ai,Bi}. 
So we can deduce the behaviour of the couphngs from our knowledge of the scahng 
of ujj. This is derived in the Appendix. We find that the couphngs are regular as 
a — >• 0, except if k < [n/2] in which case, 

n _ _ _ 

Vf ^ (A'.a-"/-)-'5:Res.. ™ . (3.18) 

where the {uji} are the one-forms on E_. Therefore, in the double-scaling limit 
proposed in ( ^.12| ), we find that these interactions remain finite ~ A^^, while the 
other 3-point vertices — * 0. This is yet further evidence of the decoupling of the 
Hilbert space into two decoupled sectors where the interactions in the 7i_ sector 
remain finite in the double-scaling limit while those in 7i_|_ go to zero. Notice, also 
that these interactions of the 7i_ sector depend universally on S_. 

The final F-term quantity that we consider is the Hessian matrix for the glueball 
superfields 



Using ( |3.17| ) we find 

< = E N.^^^ = -L V Res, , (3.20) 



T-N^t'-l^i. (3.21) 



where we have defined the 1-form T 

It is known that T can be can be written simply in terms of the on-shell Seiberg- 
Witten curve [39]: 

T = d\og{PN + ysw) . (3.22) 

In the limit a — > 0, we can take the near-critical expressions for ysw in (|2.28| ) and 
for Pn{x) = 2A^ to get the behaviour 

T — > A-'-"/2/j^,(o)_^!_a"/2^y/^(^) _ jVa"/2 . (3.23) 

We also need 



dy a'"+"/2rf(Z„(x)-\/5(5)) ~ a"^+"/2 . (3.24) 
The scaling of the holomorphic differentials is determined in the Appendix. 

Counting the powers of N and a, we find that for any j and k, HjjJ goes like an 
inverse power of a and hence diverges in the double- scaling limit (the powers of 
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cancel). This, however, presents us with a puzzle. In the case without double points 
described in [4], the Hessian was shown to vanish for the Ti- sector, i.e. j, k < [n/2]. 
Let us see how this is compatible with the scaling we have just seen. In the case, 
J,k<[n/2], 



^g)^,-(n..i)^Res, 



where 



1=1 



UJj 



d\ BJx) UjUk 



dx d[ Zm(x) \ BJx 



Lj (x) 



BJi) 



■.dx. 



(3.25) 



(3.26) 



and Lj{x) is a polynomial of degree [n/2] — 1. Note that the differential UjUk/dx 
has simple poles at x = 6/ on the curve E_: 



UJj UJk 



dx Bn{x) 
but has no pole at x = oo. For example for n odd, we find 

UJj UJk dx 



dx 



x'^ 



(3.27) 



(3.28) 



This means that in the case with no double points, m = 0, the Hessian matrix 
elements ( p.25 ) vanish identically: 



1=1 



UJj UJk 

. dx 







(3.29) 



because the sum of all residues of a meromorphic differential on the compact near- 
critical curve S_ is identically zero. This is precisely the result found in [4]. On the 
other hand, if m > 0, the Hessian matrix element will not vanish in general, because 
the differential on the right-hand side of ( p.25| ) has extra simple poles at the roots of 



(3.30) 



This result is very significant because it highlights an important difference be- 
tween the case with and without double points. Even though we do not have control 
over the kinetic terms of the glueball states, we take this behaviour of the Hessian 
matrix to signal that, with double points, the masses of the glueball fields become 
very large in the double-scaling limit. This is to be contrasted with the case without 
double points studied in [4], where the appearance of the [n/2] massless glueballs was 
interpreted as evidence that supersymmetry is enhanced to A/" = 2 in the double- 
scaling limit. 
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4. The Double- Scaling Limit of the Free Energy 



In this section, we will consider the behaviour of the free energy in the limit a 
0. The most powerful methods for calculating the Fg involves using orthogonal 
polynomials which we review in Section 4.2; however, as we shall see, these techniques 
appear only to be successful for the cases with one or two branch points and any 
additional number of double points. The only known way to calculate the Fg in 
general involves analyzing the loop equations and in particular using the algorithms 
recently developed in [40,41]. We use these techniques to prove our scaling ansatz 
( |2.12| ) in Section 4.3. Before this, it is useful to consider a situation which is exactly 
solvable. This is the case when there are two branch points and no double points: 
n = 2 and m = 0. 



4.1 Simple example: two branch points 

In this section, we consider in some detail the case when the critical point involves the 
collision of two branch points. Most of the formulae that we use are taken from [42]. 
With two branch points and no double points, we can engineer as in Section |^ with 
a curve of the form 

y^ = (x^-a^)(x^-b^) (4.1) 
in the limit a — > 0. In order to be on-shell, we require 

6^ = + 4A2 (4.2) 

because then 

which is the Seiberg-Witten curve of an SU (2) theory. The potential required is 

W{x) = Ne{-x^ - (a^ + 2A'^)x) . (4.4) 
3 

As a 0, by hypothesis, there is a light dibaryon with mass controlled by 

S = N f V(a;2-a2)(x2-a2-2A2)cia; 

-^-^ (4.5) 



-NVa^+2A^{2{a^ + A^)E{ay{a^ + 2A^)) + 2A^K{ay{a^ + 2A^))) . 
3 



At leading order for small a 



S = -^NAa^ + ■■■ (4.6) 
v2 
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Notice that the double-scahng hmit involves N ^ oo with Na'^ fixed, in agreement 
with ( p.llj ). One way to "see" the dibaryon is to calculate the low energy couplings 
of the two f/(l)'s. The coupling constant matrix has the form 

where r is the period matrix of the curve which can be written explicitly in terms of 
elliptic functions as [42] 

_ K{ay{a' + 

E(2AV(a2 + 2A2)) • ^ ' 

In the limit a ^ 0, the leading-order behaviour is 

(4.9) 



log(a/a2) 



for some constant a. This kind of logarithmic running of the coupling is characteristic 
of a one-loop effect of a particle of mass ~ S which is magnetically charged under 
the two f/(l)'s. 

However, in this case we can also extract the behaviour of the higher genus Fg 
as a ^ 0. In the present case, near the critical point we can approximate the curve 
by 

^ AA\x^ - a^) . (4.10) 

Now we can appeal to universality in the double- scaling limit. The reduced curve 
above is precisely the one encountered in the matrix model with a Gaussian potential 

W{x) = Nekx^ . (4.11) 

So the free energy of our original matrix model in the double-scaling limit should be 
equal to the free energy of the Gaussian model with the 't Hooft coupling S = ggN 
of the latter identified with S above. The Gaussian matrix model can be solved 
exactly [25] yielding 

(4.12) 

Hence, by appealing to universality, we deduce that in the double-scaling limit of the 
original theory as a — 

which verifies in a simple example the scaling hypothesis (|2.11|) . In particular, notice 
that the logarithm in Fq is directly attributed to a one-loop renormalization of the 
coupling due to the light state. 
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In [43], the methods of [40,41,64] are used to evaluate the genus one and two 
terms of the matrix model free energy for the class of double-scaling limits considered 
in [4]. The simplest case corresponds to the above conifold singularity with two 
branch points colliding and no double points, n = 2,m = 0, and the results match 
the genus zero, one and two terms in ( [4.12| ). This is precisely due to the fact that the 
near-critical spectral curve reduces to the curve of a Gaussian matrix model, which 
is simply a Riemann sphere. 



4.2 Orthogonal polynomials. 



Motivated by the study of random surfaces and 2-d quantum gravity, the double 
scaling limits of multi-cut matrix models were investigated in the early 90 's using 
orthogonal polynomials. This technique was originally developed in [44] as an alter- 
native way to calculate the matrix model free energy 

exp (F) = / ci<l exp - ^;^Tr W{^)^ = dxi A\x) e "f^"' ^ (4.14) 

^ ^ i=i 

where A(x) = Hj^jl^j ~ ^j) is the Vandermonde determinant, and its main features 
can be found in the reviews [45,46]. The basic idea is to calculate the integral 
in ( |4.14| ) by introducing the set of poljTiomials orthogonal with respect to the 



measure J dx e ^"^^'^ -P„(a;) Pmix) = hnSmn, where P„(x) is a polynomial of degree 
n normalized such that Pn{x) = a;" + ■ ■ ■ . Orthogonal polynomials satisfy recursion 
relations of the form xPn{x) = Pn+i{x) + crnP„(x) + r„P„_i(x), where r„ and cr^ 
are x-independent purely numerical coefficients, and cxn = if the potential is even. 
They can be determined by solving the non-linear recursive equations 

ngs hn-i = j e -^''^^^^ W\x) P„(a;) P„_i(x) , (4.15) 
= y" e-3="'^(")iy'(x)P„(x)P„(x). (4.16) 

In the large- iV limit, the re-scaled index n/N becomes a continuous variable ^ G [0,1], 
and the orthogonal polynomial method becomes useful provided that we have an ap- 
propriate ansatz for the large-ra behaviour of r„ and (T„. In the double scaling limit, 
the method provides important information about the resulting models. The reason 
is that the free energy turns out to be described in terms of particular solutions 
of specific classical integrable hierarchies. The relevant solutions are singled out by 
additional equations known as "string equations". This important result permits 
the calculation of the complete perturbative (topological) expansion of the free en- 
ergy in terms of the double scaled parameters, and it also provides non-perturbative 
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information about the resuhing models. Moreover, the relationship with integrable 
hierarchies makes it possible to write flow equations that interpolate between differ- 
ent models. 

The simplest large n ansatz is r„ — > r\^]S) and cr„ a{C,; S), such that both 
coefficients become smooth functions of ^ and the 't Hooft coupling S = QsN. It 
corresponds to saddle point configurations whose large- resolvent has just one cut 
[44,47]. In this case, to describe all the resulting double scaling limits, it is enough to 
restrict ourselves to the case of Hermitian matrix models with even potentials. Then, 
the topological large- iV expansion of the free energy exhibits a singular behaviour of 
the form 

F = J2 ^9(3) ' PgiS) -(S- Scf^^^^^^-'^ , m > 1 , (4.17) 

up to a few regular terms in Fq. In this approach it is customary to keep the potential 
W fixed and consider just the dependence on S; however, we could equivalently fix 
the value of S = Sc and look at the critical behaviour as a function of the coupling 
constants in W.^ Adopting the first point of view, both Sc and the integer m depend 
on W. In the saddle point method, the critical behaviour ( [4.17| ) corresponds to 
the case when m double points come together with one of the branch points of the 
resolvent, which can be realized with potentials of degree larger or equal 2(m + 1). It 
leads to the double scaling limit N 00, S ^ Sc, with N{S — Sc) 2(m+i) = finite. 
Using orthogonal polynomials, this limit can be performed with a scaling ansatz of 
the form 

r{^;S) = rc(l + -!^ + ---) , a(e;5) = 0, (4.18) 

/ ^ _2m+2\ 

together with ngg = ^S = Sc [1 ~ z N 2m+3 j . It is worth noticing that, in the 
saddle point approximation, this ansatz corresponds to a resolvent with two branch 
points symmetrically located at [44,47] 

x = ±Xb = ±Vr(l; S) ~ ±{xc + 0(iV~ra)) . (4.19) 

Then, 

a = Xb - Xc ^ {1 - S/Sc) ^ = , (4.20) 

and the double scaling limit prescription can be written as Na^~^2 = const., which 
is in agreement with eq. (2.12), and ensures that the free energy scales as in ( p.llj ) 
with n = 1. The resulting models are well known [45]. The free energy is given by 
(PF/dz"^ = u, where u = u{z) is a particular solution of the KdV hierarchy that 

^This is exactly what it done in the Dijkgraaf-Vafa matrix model, where S is fixed by the 
requirement to extremize the glueball superpotential. 
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satisfies a string equation depending on m. In the simplest case, m = 1, the string 
equation is 2; = m^ + m"/3, which is known as Painleve I; the resulting model describes 
pure 2-d gravity. For m > 1 the model corresponds to the (2m + 1,2) CFT minimal 
model coupled to 2-d gravity. It can be shown that the m-th string equation is 
associated to the m + 1-th flow of the KdV hierarchy. Although in the matrix model 
approach m > 1, the integrable hierarchy formulation allows one to construct an 
additional model associated to the 1^* flow of the hierarchy. Actually it corresponds 
to a topological theory that underlies all the other double scaled models, and it was 
identified with 2-d topological gravity [48]. 

Saddle point configurations with two cuts are recovered with a "period-two" 
ansatz; i.e., by assuming that r2n, r2n+i, ^2x1 and a^n+x approach different smooth 
functions in the large-n limit [49]. However, the resulting cuts always open up around 
minima of the potential W[x) [47,50] and it is not known how to describe saddle 
point configurations with cuts open around the maxima of W[x) using orthogonal 
polynomials. This includes, for instance, the two-cut phase of the matrix model 
specified by 1^ = \gj,x^ + g\X, which was extensively discussed in [42,51], as well as 
in section 4.1. Actually, this limitation poses a serious restriction for the use of the 
orthogonal polynomial method in the context of (|1.1CI| ). In any case, it is interesting 



to review the main results achieved by considering the class of two-cut matrix models 
that can be studied using this method. 

Two-cut Hermitian matrix models with even potentials were investigated in [52]. 
They exhibit critical behaviours corresponding to cases where two branch points come 
together with an odd number of double points, and lead to the same double scaling 
limits that were previously found in the study of unitary matrix models [53], which are 
described in terms of the modified KdV hierarchy. The interpretation of the resulting 
models in terms of super CFTs coupled to 2-d supergravity was suggested in the early 
90's [52,54], and it has been recently clarified and supported by Klebanov, Maldacena 
and Seiberg [55]. However, these are not the only two-cut models that have been 
studied using orthogonal polynomials. Two-cut Hermitian matrix models with both 
even and odd terms in the potential were considered in [54,56,57]. Remarkably, 
the corresponding double scaling limits are described in terms of several integrable 
hierarchies associated to s/(2, C). Following [57], the results can be summarized as 
follows. First, consider the case of real potentials. In the saddle point approximation, 
their critical behaviours correspond to configurations where an odd number of double 
points come together with two branch points [54,57]. Their double scaling limits are 
described in terms of solutions of the non-linear Schrodinger (NLS) hierarchy singled 
out by string equations associated only to the "even" flows of the hierarchy. 

However, as pointed out in [57], the most general critical behaviour, where an 
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arbitrary number of double points collide with two branch points, is obtained by 
considering complex potentials. Namely, the potential has to be taken such that 
W*{x) = W{—x), which makes the coefficients (T„ pure imaginary. This gives rise 
to double scaling limits described in terms of solutions of the Zakharov-Shabat (ZS) 
hierarchy specified by string equations associated to the complete set of flows, but 
the ffist one. The resulting models are described by two real functions f{z) and g{z). 
They enter the ansatz for the orthogonal polynomial coefficients, which is of the form 

l + i-ir4^ + ---) , an = tajb+{~ir4^ + ---) , (4.21) 

J\[ m+2 / \ ]\[m+2 / 

/ ^ m + 1 \ 

together with = Sc [1 — z N '"+2 j . Here, b is an arbitrary number, and m is the 
number of double points. In the saddle point approximation, this ansatz gives rise 
to a resolvent where the two colliding branch points are located at [47, 50] 



^ (T2n + o-2„+i ± y (cr2„, - cr2„+i) - 4(^7^ - ^r2„+i) 

ibac±0{N'^) . (4.22) 

Thus, in this case, 

a = x^ -x; {l- S/Sc) ^ = , (4.23) 



and the double scaling limit prescription can be written as Na^^"^ = const., which 
is again in agreement with eq. (2.12) for n = 2. The free energy is given by 
cPF/(fz = [f^{z) — g'^{z)) /4. For m = 1 the string equations are 

zf + /(/ - n + 2f" = , zg + g{g^ - f) + 2g" = , (4.24) 

and for m = 2 they read 

2zf + 3g\g' - f) + 2g"' = , 2zg + Sf'ig' - f) + 2/"' = . (4.25) 

They exhibit that the string equations corresponding to odd m admit the reduction 
g = and the analytical continuation g ig, which lead to the mKdV and the 
NLS hierarchies, respectively. Other reductions, as well as the relationship between 
the partition function and the tau-functions of the hierarchies, have been discussed 
in [57]. It is straightforward to check that the free energy scales exactly as in ( p.llj ) 
with n = 2. 



^°This is equivalent to consider anti-Hermitian matrix models, whose eigenvalues are purely 
imaginary, with real potentials. 
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An interesting observation about the structure of these two-cut matrix models 
was made in [54]. There, it was shown that the formulation in terms of the ZS 
hierarchy allows one to define a model associated to the 1*^* flow which is of topological 
nature. This is similar to what happens in the one-cut case, where the 1®* flow of 
the KdV hierarchy leads to 2-d topological gravity. By comparison, it is natural to 
expect that the topological model of [54] underhes the structure of all the ZS-related 
two-cut matrix models, although the proper interpretation of this topological phase 
is still unclear. 

Another important feature pointed out in [50] and, mostly, in [55] is the fact 
that the string equations of these models admit an integration constant. For the 
models associated to the ZS hierarchy, this is a consequence of the invariance under 
the S0{1, 1) transformations f ± g ^ e ^^{f ± g). This constant shows up in the 
perturbative expansion of the free energy, although the leading order term remains 
independent of it. The only exception is the topological model of [54], where the free 
energy is proportional to a non- vanishing integration constant. However, in [54] the 
presence of this constant in the non-topological models was missed, mainly because 
it was assumed that the potential is real and, consequently, the colliding double 
points and branch points are on the real axis. In [55], double points were allowed to 
be complex, the integration constant was understood in terms of the period of the 
hyper-elliptic y = y{x) curve of the matrix model around the cut joining the two 
colliding branch points, and its physical interpretation was clarified. When g — > ig, 
and the ZS hierarchy changes to the NLS hierarchy, the constant becomes complex 
too, and it can be related to the difference between the number of eigenvalues sitting 
on the two cuts [50,55]. 

Motivated by the correspondence between the period-2 large-n ansatz and two- 
cut saddle point configurations, a study of multi-cut configurations by considering 
more general period-g large-n behaviours; i.e., riq+i — > ri{^; S) and cx/g+i •S'), 
for i = 0, 1, . . . , g — 1, was initiated. However, for g > 3 this sort of ansatz is 
only able to reproduce multi-cut saddle point configurations for very special types 
of potentials. For instance, ior q — 3 and even potentials of degree six with three 
minima, the correspondence turns out to works only in the particular case when 
all the minima are degenerate [47]. In other words, the appropriate large-n ansatz 
corresponding to generic saddle point configurations with three or more cuts in not 
known, which prevents the use of the orthogonal polynomial method to investigate 
the properties of these configurations and their double scaling limits. Several aspects 
of this problem were discussed in [58] . 

An interesting particular set of multi-cut matrix models partially studied using 
orthogonal polynomials are provided by the, so-called, "orbifold (Hermitian) matrix 
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models" proposed in [52]. They give rise to saddle-point configurations with 2k 
cuts symmetrically radiating from the origin. The critical behaviours considered in 
that paper correspond to cases where the end points of all the cuts simultaneously 
collide together with 2km — 1 double points, for m > 1. The resulting double- 
scale models are expected to exhibit a symmetry. The simplest potential that 
gives rise to such configuration is W{x) = — x'^ for a matrix model defined over 
matrices $ such that is hermitian, for which there are four branch points and 
three double points in the critical region. This model was studied in [52] by means 
of a very symmetrical period-4 large-n ansatz for the coefficients r„. However, both 
its physical interpretation and the associated integrable hierarchy remain mysterious 
still. Nevertheless, it can verified from the resulting generalized string equations, 
that the free energy scales as in (|2.11|) with n = 2k and m = 2km — 1. 



Apart from the aforementioned models, the fact that the orthogonal polyno- 
mial method seems to work for one or two branch cuts may be related to the fact 
that resulting near-critical curves are both of genus zero. In fact, this is related 
to the fact that the string equations in both cases can be formulated in terms of 
the Virasoro algebra of a chiral boson on the complex plane. For the case of one 
branch point there is a branch cut extending to infinity and the boson is twisted. 



4.3 The loop equations 



In this section, consider a radically different approach for investigating the free energy 
which involves a recursion relation for the Fg. Before we come to the loop equation 
itself we need to define the p-loop correlator, or p-point loop function, as 



W{xi, ...,Xp)= NP-\ Tr ^ ■ ■ • Tr 



Xi — $ Xp — ^ ' conn 



(4.26) 



It has the following genus expansion 

W{xi, ...,Xp) = J2 TvT^^'Ha^i, ■■■,^p)- (4.27) 

In [40], Eynard found a solution to the matrix model loop equations that allows to 
write down an expression for these multi-loop correlators at any given genus in terms 
of a special set of Feynman diagrams. The various quantities involved depend only 
on the spectral curve of the matrix model and in particular one needs to evaluate 
residues of certain differentials at the branch points of the spectral curve. 
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This algorithm and its extension to calculate higher genus terms of the matrix 
model free energy [41] represent major progress in the solution of the matrix model 
via loop equations [51,59-64]. This is particularly important because, as we have 
seen in the last section, the orthogonal polynomial approach can only applied in 
very special cases. A nice feature of the loop equation approach is that they show 
directly how the information is encoded in the spectral curve. In particular, we will 
be able to make some precise statements on the double-scaling limits of higher genus 
quantities simply by studying the double-scaling limit of the spectral curve and its 
various differentials. 



Given the matrix model spectral curve for an s-cut solution in the form ( p.4|) the 
genus zero 2-loop function is given by 



W{Xi,X2) = 



cr'(si) A{xi,X2) 



(4.28) 



4:{Xi- X2)^/oiJ^^/oix2) 4a/o-(Xi) a/(t(x2) 

where A is a symmetric polynomial given by 

Aix,,X2) = f:^^^^^^^^, (4.29) 

with 

^-2 s-1 r dx 1 
A(X2) = 5^A^4 = -E^^(^2) / T (4.30) 

and s is the number of cuts. The polynomials Lj{x) are related to the holomorphic 
1-forms and defined in the Appendix. 



The genus zero 2-loop function for coincident arguments is 

TT^/ \ 1- TT^/ \ cr"(xi) cr'(xi)^ A(xi,Xi) 
Wix., X.) ^ Jmr^ (x„ x.) ^ + + 

1 o'l Ux) 



(4.31) 



16(x - ai)2 16a--(x - ai) 4(x - ai) 
The other important object is the differential 

dS2i-liXi, X2) =dS2i{Xi,X2) 



(r(X2 



A/cr(xi) \X1 - X2 a/(t(x2) 



^i{x\) _ ^Q(^2;2)Lj(a;i) j (ixi , ^ ^ 
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where i = 1, . . . , s and 



A crucial aspect of the one-form ( [4.32| ) is that it is analytic in X2 in the limit X2 

(J2i-i or (T2i [40] 



dSi{xi,X2) 

lim 



(7(xi) (2:1 -X2 ^^^^^^^^ Ia, ^/^{x~X2) ) 



(4.34) 

The subtlety is that in the definition of ([4.33|) , the point X2 is taken to be outside 
the loop surrounding the j-th cut, whereas in ([4.34|) , X2 is inside the contour. Note 
also that 

A{x^,X2) = -X^ (j2Lj{^2)C,{a,)] (4.35) 

1=1 \j=l / 

and in particular 

A{xi,a,) = C,{xi)a'{ai) . (4.36) 



The expression of W^^\xi, . . . be found by evaluating a series of Feyn- 

man diagrams of a cubic field theory on the spectral curve [40]. To this end, define 
the set Tp^^ of all possible graphs with n external legs and with g loops. They can 
be described as follows: draw all rooted skeleton trees ( trees whose vertices have 
valence 1, 2 or 3 ), with p + 2g — 2 edges. Draw arrows on the edges from the root 
towards the leaves. Then draw in all possible ways p — 1 external legs and g inner 
edges with the constraint that all the vertices of the whole graph have valence three, 
namely that are always three and only three edges emanating from any given vertex. 
Each such graph will also have some symmetry factor [40]. 

Each diagram in then weighted in the following way. To each arrowed edge 
that is part of the skeleton tree going from a vertex labelled by xi to a vertex 
labelled by X2 associate the differential dS{xi,X2) (H^.32|) . To each non-arrowed edge 
associate a genus zero 2-loop differential G{xi,X2) = W{xi,X2) dxidx2 and to each 
internal vertex labelled by xi associate the factor {2eNy{xi)dxi)~^ . For any given 
tree T e Tp^\ with root xi and leaves Xj,j = 2, . . . ,p and with p + 2g — 2 vertices 
labelled hj x'^,v = 1, . . . ,p + 2g — 2, so that its inner edges are of the form vi V2 
and its outer edges are of the form v —>■ j, we define the weight of the graph as 



- 26 - 



follows 

P+2g-2 2s 

= j^j^, n E ^5,^ n ^S:i<. <■) 

v=l iv=l \ u/ u jjjj-igj. g(jggs tj^u, 

inner non-arrowed edges v'^w' outer edges t)— >j 

(4.37) 

In order to find an expression for Fg, g > 1, one should consider the same graphs 
relevant for W^^\xi) and do then the following [41]: 

(i) Eliminate the first arrowed edge of the skeleton tree. Labelling the first vertex 
by Xi and the second vertex by X2, this amounts to dropping the factor dS{xi, X2). 



(ii) The factor {2eNy{x2)dx2) ^ has to be dropped and replaced by 

f^^ y(s)ds 

; . (4.38) 

y{x2)dx2 

Note that when evaluating the final residues at X2 = (Ji, one needs to expand the 
above integral by setting go = [41]- It is also understood that the evaluation of 
the residues starts from the outer branches and proceeds towards the root. This 
procedure does not apply for the genus one free energy whose expression has in any 
case been found via the loop equations in [51,63,64]. 



We will consider the a ^ limit of each element in ( [4.3?| ). Using the results in 



the Appendix for the scaling of L, it is straightforward to argue that for a branch 
point bi in the critical region 

dSi{xi,X2) — > dSi{xi,X2) 



2_^Cj[x2)Lj[xi) dxi , 



where dSi is the analogous differential on E_ and Lj{x) — * a"'^'^~^Lj(x) for j < [n/2]. 
Conversely, the differentials dSi{xi,X2) where i labels a branch point of the spectral 
curve that remains outside of the critical region give a vanishing contribution in the 
double-scaling limit. Likewise using equations ( |4.28| ), ( [4.29| ), ( ^4.301 ) and ( |4.35| ) we 
have 

G{xi,X2) = W{xi,X2) dxidx2 — > G{xi,X2) = W{xi,X2) dxidx2 , (4.40) 
where G{xi,X2) is exactly the 2-point loop correlator on S_. 
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So far we have seen that the double points of the near-critical curve do not play 
a role in taking the limit of the differentials. However, this is not the case for the 
final two elements of the Feynman rules 

ydx — > v^a'"+"/2+i y_ dx (4.41) 

and 

I,^y(^)ds _^ /; y-{~s)d~s ^ 

y{x)dx y_{x)dx 
To summarize: what we have found is that all the relevant quantities reduce to the 
analogous quantities on the near-critical curve in the limit a — 0. In particular, 
being careful with the overall scaling, the genus g free energy has the limit 

Fg — ^ c^-9/s^2-^g Fg{^_) . (4.43) 
where we have emphasized that Fg(E_) depends only on S_. This is the result 



advertised in (|2.11|) and the property of universality. Similarly, the genus g p-point 



loop functions have the limit 

Wg{Xi,..., Xp) dXi--- dXp ^ C^^9-p/2^2~2g-p ^^^(^^^ . . . ^ ^Xi ■ ■ ■ dXp . (4.44) 



5. The Physics of the Double Points 



In this section, we consider in more detail the case when the near-critical curve has 
double points, since as we have already argued the physics is very different from that 
considered in [4] with only branch points. 

It will turn out that the existence of the double points has a profound effect on 
the physics of the critical region. In the case of branch points, dibaryons become 
massless when the branch points collide. We want to argue that certain states also 
become massless when a double point collides with a branch point. In particular the 
mass scale when m double points and n branch points come together is Na"^'^"'^'^'^^ . 
In order to tease out the physics of these light states let us consider the case with 
only a single branch point and one double point. The curve in the near critical region 
has the form 

y^^ = (^x - zYix ~b) . (5.1) 

This is precisely the kind of critical curve that appeared in the old matrix literature 
[45]. 

In order to motivate the issues involved, let us consider what happens as we go 
through a transition where the double point arises from the collision of two branch 
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points. This can be achieved by using a more general form of the bare superpotential 
than the ones we considered in Sections 2.2 and 2.3. In this case, on one side of 
the transition we have three branch cuts ending in the critical region and three 
light dibaryons associated to the three cycles Ci which link each pair of branch 
points. Each of the cuts is associated to an abelian gauge field U{l)i in the IR gauge 
group. The situation is illustrated in Figure The dibaryon Ci is magnetically, 




\ \ 
\ \ 
\ \ 

U(1), 

(a) (b) 

Figure 1: (a) A configuration for which S_ consists of three branch points, (b) Two of the branch 
points coalesce to form I]_ with a single branch point b and a single double point z. 

and potentially electrically, charged with respect to U{l)j x U{l)k- In fact, denoting 
the electric and magnetic charges as {gi,qi', 92,92] 93,93) we can choose conventions 
where the three dibaryons have charges 

Ci: (0,0; 1,1; -1,0) 

C2: (-1,0; 0,-1; 1,0) (5.2) 
C3: (1,0;-1,0;0,0) . 

The three dibaryons have mutually non-local charges since as cycles Ci-C2 = C2-Cs = 
C's-C'i = 1 and C1 + C2 + C3 = [17]. Now suppose two branch points come together 
corresponding to the vanishing of C3. At this point the dibaryon C3 becomes massless 
and the double points of the M = 1 and Seiberg-Witten curves coincide. This point 
in the parameter space touches a new phase where the dibaryon is condensed — the 
double points of the two curves move apart — and the U{1) subgroup of f/(l)i x f/(l)2 
defined by the elements (e*^, e~*^) is confined. This new confined phase is the one 
we are interested in. Since dibaryons Ci and C2 carry opposite electric charge under 
the confined U{1) it is natural that Ci and C2 form a bound state: a "meson". The 
bound state has minus the charges of C3; in other words, the meson is magnetically 
charged under the confined U{1). Notice, however, that the state is neutral under the 
two remaining [/(l)'s. This means that, in stark contrast to the case of dibaryons, 
we will not "see" these bound states via the logarithmic running of the coupling of 
the IR abelian gauge fields as we did for the case of a light dibaryon in Section ^ 
However, as we have seen, their effects will show up in other quantities. In particular, 
if this meson, of mass M, runs around 1-loop graphs then it contributes as M^~^^ 
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to the higher-genus free energy Fg. So the fact that Fg ~ means that since 

m = 1 and n = 1 the mass of the meson is M ~ A = Na^^'^, in this case. 

The remaining issue is to argue that the meson becomes massless when the double 
point coincides with the remaining branch point, since as a — evidently M ^ 0. 
One might have thought that the meson has a mass of order the confinement scale 
and would not become massless at the critical point. The issue is rather reminiscent 
of the pion in QCD. The pion mass is much lower than the scale of confinement 
because the pion is the would-be Goldstone boson of chiral symmetry breaking. 
However, chiral symmetry is explicitly broken by the mass of the light quarks and so 
a mass for the pion is generated but at a scale much smaller than the confinement 
scale. In the present situation, when the double point meets the branch point and 
the Seiberg-Witten double point, z = b = h, the dibaryon condensate vanishes and 
the theory is at an A/" = 2 Argyres-Douglas superconformal field theory [17]. As the 
Seiberg-Witten double point moves away (but leaving the branch point and double 
point of the A/" = 1 curve coincident h ^ z = h) the dibaryon condensate develops 
and M = 2 superconformal symmetry is spontaneously broken. As the Seiberg- 
Witten double point moves out to infinity, only an A/" = 1 superconformal symmetry 
remains. Our proposal is that the mesonic bound-state is the Goldstone mode of this 
broken symmetry and hence is exactly massless when the branch and double points 
coincide. As the double point of the M = 1 curve moves away from the branch point, 
z ^h, this state becomes massive with a mass ~ A as indicated by the behaviour of 
Fg in the double-scaling limit ( |4.43| ). 

In the general case, there will be bound states of the type described associated 
to each pair consisting of a branch point and a double point. These bound states will 
be magnetically charged with respect to the confined U{1) associated to each double 
point. In the limit we are considering with hi finite as a ^ 0, the confinement scale 
is always much greater than the masses of the bound states. 



6. Conclusion 

In this paper, we have studied the large- A^ limit of certain M = 1 theories in the 
proximity of Argyres-Douglas-type singularities. The exact analysis of the F-terms 
performed via the Dijkgraaf-Vafa matrix model correspondence shows that, close to 
these singular points, the 1/A^ expansion breaks down. This can be traced back to the 
appearance of extra light particles of baryonic and/or mesonic nature that become 
massless at the critical point. As in [4], there is a natural large- A^ double- scaling limit 
that emerges: namely if one considers approaching the singularity in conjunction with 
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the 't Hooft large- limit in such a way that the mass M of these baryonic or mesonic 
states is kept fixed. In this limit, we have argued that the Hilbert space of the theory 
splits into two mutually decoupled sectors, 7i+ and The sector ?i+ obeys the 
usual large- scaling and becomes free in the limit, whilst the sector 7Y_ keeps non- 
trivial interactions weighted by the effective string coupling (yfes ~ 1/M. In [4], it 
was possible to map this large- iV double-scaling limit to an analogous double-scaling 
limit considered in [12] in the context of the duality between four- dimensional Little 
String Theories and certain non-critical string backgrounds [15]. This led to the 
proposal that the non-trivial dynamics of the 7i_ sector admits a dual holographic 
description in terms of the four-dimensional non-critical string. This proposal passes 
a non-trivial consistency check in that the non-critical string dual exhibits M = 2 
supersymmetry while the F-term effective action was shown to be consistent with 
an enhancement from M = 1 to M = 2 supersymmetry in the 7i_ sector. In this 
paper, we have shown that for more general singularities the sector 1-L_ has only 
N = 1 supersymmetry. In particular, the effects of the glueball superpotential do 
not vanish in the limit. Thus, even though the field theory analysis clearly shows 
that the double-scaling limit still selects a particular sector of the theory which has 
non-trivial dynamics, we do not identify the would-be non-critical string dual as 
in [4]. For the type of singularities studied in the old matrix model to describe two- 
dimensional gravity coupled to c < 1 conformal field theories, where the near-critical 
curve has only one branch point along with any arbitrary number of double points, 
the double-scaling limit describes a sector of a gauge theory with a mass gap and 
light meson-like composite states which we argued are the approximate Goldstone 
bosons of superconformal invar iance. These cases are special in that there are no 
massless abelian fields in the 7i_ sector; in fact, there is a mass gap. 
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Appendix A: Some Formulae 



In this appendix, we consider the double-scaling limit of various quantities de- 
fined on the curve S ( p.4| ). This is most conveniently done in the basis {Ai^Bi} of 
1-cycles described in Section ^. In particular, for i < [n/2\ these are cycles on the 
near-critical curve S_ in the double-scaling limit. 
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The key quantities that we will need are the periods 

Mij = (f ^^^==dx , N,, = I ^^^==dx . (A.l) 
Jb, ^/(T{x) J A, V^(^) 

First of all, let us focus on N^j where j < [n/2\, but i arbitrary. By a simple scaling 
argument, as a — 0, 

iV., = f'' dx a--'- dx = a-«/^ fiP(h,) , (A.2) 

for some function f^^^ of the branch points of S_. Here, 6^.^ are the two branch 
points enclosed by the cycle Aj. A similar argument shows that Mij scales in the 
same way: 

M,, ^ a'-"/' flf\k) . (A.3) 

So both Nij and M^j, for i,j,< [n/2], diverge in the limit a ^ 0. On the contrary, 
by using a similar argument, it is not difficult to see that, for j > [n/2], Nij and M^j 
are analytic as a — >■ since the integrals are over non-vanishing cycles. 

In summary, in the limit a — 0, the matrices and M will have the following 
block structure 

where by — or + we denote indices in the ranges {1, . . . , [n/2]} and {[n/2] + l, . . . , s — 

1} respectively. In ( [A. 41 ), A^ and are divergent while the remaining quantities 

are finite as a — >• 0. 

We also need the inverse L = N"'^. In the text, we use the polynomials Lj{x) = 
Ylk=i^jk^^~^^ which enter the expression of the holomorphic 1-forms associated to 
our basis of 1-cycles, 

(t ujj = 6ij . (A. 5) 

These 1-forms are equal to 

uj.{x) = -^^ikL dx = ^kJiLjkx'-' ^ r ^^^^^ ^ 

where j = 1, . . . , s — 1. From the behaviour of A^ in the limit a — > 0, we have 

L = AT-i -> ( - ^^^^ -1 ^ , M= -N-J_ N^l (n^I) " . (A.7) 
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Since is singular we see that L is block diagonal in the limit a ^ 0. This is just 

an expression of the fact that the curve factorizes S — > S_ U S+ as a — >■ 0. In this 
limit, using the scaling of elements of Lj^, we find, for j < [n/2]. 



^["/2l^ f(N)\-l ~.k-l 



s-^vniz\, f(N)\- 

ujj — > ax = ojj . (A-oj 



'B{5:) 

the holomorphic 1-forms of E_. While for j > [n/2]. 



^^>[V2](^++)^^ (A.9) 



are the holomorphic 1-forms of S 



+ ■ 
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